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We discuss the physics underlying an all-order resummation of logarithmic en¬ 
hanced contributions to dijet cross sections, and present preliminary results for 
the distribution in the dijet transverse energy difference in DIS. 


1. Introduction 

Clustering hadrons into jets is a very useful tool to make QCD predictions. 
Measuring jet cross section removes the theoretical uncertainty due to frag¬ 
mentation functions, and makes it possible to directly compare data with 
perturbative (PT) QCD predictions, whose degrees of freedom are partons 
and not hadrons. This comparison works extremely well when considering 
high transverse energy {Ep) jets. In this kinematical situation one is ex¬ 
pected to probe quarks and gluons at small distances. This joint theoretical 
and experimental effort has led to the measurement of the QCD coupling 
Us from jet inclusive Ep spectra and to better constrain the gluon density 
in the proton^. 

There are actual limitations in the use of perturbation theory. First of 
all any perturbative series diverges. This is related to the fact that the 
observed degrees of freedom are hadrons and not partons, and results in 
an ambiguity in PT predictions that fortunately is suppressed by inverse 
powers of the hard scale of the process. The second limitation has to do with 
the fact that the coefficients in the PT expansion can be logarithmically 
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enhanced due to incomplete real-virtual cancellations in the infrared (IR). 
This has the consequence that in particular phase space regions one observes 
a breakdown of the PT expansion, which can be cured only by performing 
an all-order resummation of such logarithms. 

It is therefore mandatory, in order to have predictions valid in the whole 
of the phase space, to complement any fixed order calculation with the all¬ 
order resummation of logarithmic enhanced contributions, and, if possible, 
to remove the ambiguity in the PT expansion with a power-suppressed non- 
perturbative (NP) correction, which, unless it can be computed from the 
lattice, has to be taken as a phenomenological input^. 


2. The observable 


The observable we consider is the dijet rate, i.e. the fraction of events with 
at least two jets. To select hard jets we put a cut on the transverse energy of 
the two highest Et jets, requiring Eti > Et 2 > E^- This particular choice 
is referred to as “symmetric Et cuts”. It was noted by Klasen and Kramer^ 
that symmetric Et cuts produce IR instabilities in next-to-leading order 
(NLO) QCD predictions. It was later Frixione and Ridolfi"^ who proposed 
to perform an asymmetry study by considering 

(t(A) = a{Eti > Em + A; Et2 > Em) ■ (1) 


They computed cr(A) in photoproduction at NLO, and obtained that while 
(t(0), the total dijet rate with symmetric Et cuts, is finite, the slope of the 
curve cr'(A) = da/dA diverges for A = 0. This behaviour is not present 
at all in the data, as one can see from the asymmetry study performed by 
ZEUS in DIS®, and reported in fig. 1. There the data is plotted against NLO 
QCD predictions obtained from the numerical program DISENT®. Again, 
the slope of the NLO curve diverges for A = 0, while the data decreases 
smoothly with increasing A. This can be easily understood by noting that 
cr'(A) is just (minus) the differential distribution of the highest Et jet: 


a'(A) 


da 


dExi 




( 2 ) 


Since the latter quantity is a physical cross section, the slope cr'(A) has to 
be negative. This corresponds to what is seen in the data, and also implies 
that any turnover in the NLO curves is unphysical. 

To better understand the origin of the divergence in a'(A), we consider 
the value of the jet transverse energy difference Eti — Et 2 - In the specific 
case of dijet production in DIS, with an incoming quark p, at the Born level, 
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ZEUS 



Figure 1. The dijet rate C7(A) as a function of Eti in DIS. The dots are the data, while 
the lines are NLO predictions. 

a dijet event consists only of two outgoing hard partons pi and p 2 , in this 
case a quark and a gluon. Their transverse momenta pti and pt 2 are back- 
to-back, so that, defining Eti = \pti\, we have Eti = Et 2 - After emission 
of a soft gluon k not clustered with any of the two outgoing partons, we 
obtain from transverse momentum conservation Eti —Et 2 — \kx\^ where kx 
is the component of the gluon kt parallel to pti- When considering cr'(A), 
Eti is forced to lie on the line Eti = Em + A, but from kinematics Eti = 
Et 2 + \kx\- These two lines should then intersect somewhere in the allowed 
phase space Eti > Et 2 > Em, and this is possible only for < A, 
as can be seen from fig. 2. Squeezing soft radiation is known to give rise 
to large logarithms in fixed order calculations. Indeed, if we consider the 
emission of a soft gluon collinear to the incoming quark together with the 
corresponding virtual correction, we obtain 

ct'(A) = CTo(A) X (^1 - > (3) 

so that the slope cr'(A) diverges for A ^ 0. A quick solution of the problem 
can be achieved by choosing kinematical cuts such that the Born slope 
CTq(A) ^ 0 for A ^ 0. This can be achieved for instance by imposing 
asymmetric Et cuts^, or choosing a suitable range for the dijet invariant 
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Figure 2. The phase space for the dijet rate in the Eti — Et 2 plane. 

mass and rapidities^. However, the actual solution relies in the all-order 
resummation of soft gluon effects, as will be explained in the next section. 


3. Resummation 

Considering an arbitrary number of partons one finds Eti—Et 2 — \ hx \, 
where the primed sum runs on all partons not clustered with the jets. 
Introducing S{kt), the probability that = kt, we have cr'(A) = 

crQ(A)bb(A), where the K-factor bb(A) is in general given by 

H^(A) = J £kt S{kt) 0(A -\kx\)=‘^ y y sin(6A)E(6). (4) 

For emissions soft and collinear to the incoming quark, no secondary par- 
tons are clustered with pi or p 2 , so that an all-order resummation of such 
contributions yields 11(6) = exp[—In^ 6]. From eq. (4) we have then 
W (A) ~ A for small A, which roughly corresponds to the behaviour seen 
in the data. 

In general, 11(6) = exp[Lgi(asL) -I- g 2 {asL) + asg^io-sL).. .], with L = 
In 6, and we aim at next-to-leading logarithmic (NLL) accuracy, i.e. the 
knowledge of gi and g 2 - In order to perform an all-order resummation, one 
need to reorganise the PT series and collect all contributions to S(6) up to 
a given accuracy. First of all one can see that at NLL accuracy, E(6) can be 
interpreted as the probability that \ Y^[kix \ < jh. Furthermore, given 
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a variable V, a generic function of all final state momenta, multi-parton 
matrix elements and phase space constraints can be approximated by using 
the following general properties: 

(1) all leading logarithms (LL) originate from configurations of soft and 
collinear gluons for which the hardest emission ki dominates, that 
is V(ki,...,kn) ~ V(ki); 

(2) hard-collinear and soft-large-angle emissions give a NLL contribu¬ 
tion, and again one can approximate V"(fci,..., ~ V(ki); 

(3) the fact that V{ki,...,kn) ^ y{ki) needs only to be taken into 
account for soft and collinear emissions; these “multiple emission 
effects” are NLL, and can be treated either analytically with a 
suitable integral transform, or numerically with Monte Carlo (MC) 
techniques; 

(4) secondary splittings are accounted for by taking the QCD coupling 
in the physical CMW scheme, at a scale of the order of the transverse 
momentum of the parent parton, i.e. as as{kt). 

The strongest implication of these four statements is that at NLL accuracy 
soft and/or collinear (SC) partons can be considered as emitted indepen¬ 
dently from the hard parton antenna, in spite of the fact that QCD is a 
non-abelian gauge theory. This crucial simplification is known as the “in¬ 
dependent emission” approximation, and leads to exponentiation of leading 
logarithms, and factorisation of NL logarithms. If we assume the validity 
of independent emission, we obtain 

S(&) = ® , R{b) = J[dk] w{k) 0 (|fc,| - e-^^/b) , (5) 

where the ’radiator’ R{b) is (minus) the contribution to S(6) from a single 
SC gluon k emitted with probability w{k). The ratio of parton densities 
/(l/6)//((5) is a general feature of resummations involving incoming par- 
tons, and is due to the fact that the scale of the incoming parton density 
is set by the upper limit of observed parton transverse momenta, which in 
this case is of order 1 /b. 

Unfortunately, independent emission approximation does not hold for 
all variables, but only for those who are recursively infrared and collinear 
(rlRC) safe and (continuously) global®. The lack of either of these two 
properties causes one (or more) of the above statements not to be true. 

For instance, for the three-jet rate in the JADE algorithm, which is 
global but not rlRC safe, statement 1 is false, since there are LL contribu- 
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tions originated by multiple emissions, with the consequence that leading 
logarithms do not exponentiate®. 

Globalness means simply that V is sensitive to emissions in the whole of 
the phase space. This is not the case for Eti — Et 2 , since only partons that 
are not clustered with the two hard jets contribute to the dijet transverse 
energy difference. The first consequence of non-globalness is that secondary 
splittings cannot be accounted for by simply setting the proper scale for as, 
because new NLL contributions, called “non-global logs”^®, arise when a 
cascade of energy ordered partons outside the measure region emits a softer 
gluon inside. In this case phase space boundaries have to be taken into 
account exactly for all emissions, so that one has to rely on MC methods. 
Moreover, due to the complicated colour structure of multi-gluon matrix 
elements, non-global logs at the moment can be resummed only in the large 
Nc limit. 

In our case, if particles are clustered into jets with a cone algorithm, 
after having specified a procedure to deal with overlapping cones, one has 
that to NLL accuracy a jet consists of all particles that flow into a cone 
around pi or p 2 , and equation (5) gets modified as follows^^: 

m = ® X sm] , t{b) = ^ , (6) 

where S(t) represents the contribution of non-global logs, and can be com¬ 
puted with a MC procedure as a function of the evolution variable t. 

If jets are clustered with the kt algorithm^^, further complications arise. 
In this case the fact that a particle belongs to a jet depends on all other 
emitted particles. This has interesting implications for non-global logs. 
Dasgupta and Salam^^ have shown that the dominant contribution to non- 
global logs arises when the observed gluon is close to the boundary of the 
measure region. Appleby and Seymour^^ noted then that the kt algorithm 
requirement forces such gluons to be clustered with gluons outside the mea¬ 
sure region, thus reducing the magnitude of non-global logs. 

One can also ask whether for the independent emission contribution to 
S(&), the exponential form of equation (5) is still valid. This is true for cones 
(the “unclustered” case) since the phase space constraint is factorised, but 
should be checked for the kt algorithm (the “clustered” case). In order to 
answer this question we considered a much simpler variable, the transverse 
energy flow away from the jets. Given a pair of hard jets (taken back-to- 
back for simplicity), one defines a region Q away from the jet axis (which 
in e^e“ annihilation roughly coincides with the thrust axis), and the away- 
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from-jet Et flow 

ET,n = ^ ku , (7) 

ieo 

with kti the transverse momentum of the i-th jet with respect to the 
jet axis^^’^"*. The quantity T,{Q,Qq,), the probability that i?T,n < Qq, 
in the region Qq ^ Q is sensitive only to soft gluons at large angles. 
One then wishes to resum LL contributions to I]((5,Qo), whose order is 

<ln”(Q/go). 

At LL accuracy E{Q, Qq) is given by 

nQ,Qn) = En,p{t)-S{t), ( 8 ) 

where t = t{Qci) is the evolution variable defined in equation (6). Glu¬ 
ons emitted directly from the two hard partons give rise to Sn,p(t), while 
non-global logs are embodied in S{t). Appleby and Seymour^"^ assumed 
that also for the clustered case Sn,p(t) = that is the single gluon 

contribution to exponentiates. This naive expectation can be 

motivated by the fact that if emissions are assumed to be independent, 
multiple emission effects are usually relevant only for soft and collinear 
gluons (statement 3). However, since Et,q, is manifestly non-global, one 
cannot exclude multiple emission effects coming from soft gluons at large 
angles. Consider for instance two gluons fci and ^ 2 , with loi ^ uj 2 , 
inside H and ki outside. It is possible that the jet algorithm cluster ^2 
with fci, thus spoiling the exponentiation of the single gluon result^®. This 
is indeed seen when comparing the resummed expression for the differen¬ 
tial distribution a~^da/dL, with L = ln(go/g), with the NLO program 
EVENT2®. If LL are correct, the difference of the two distributions should 
go to a constant for large (negative) L. This happens for the coefficient 
of CfCac^I, indicating that both in the unclustered and the clustered case 
non-global logarithms are correctly taken into account. If one however uses 
for the clustered case the same expression for Eq^p as for the unclustered 
case, one finds a discrepancy in the coefficient of Cpol, the plot of fig. 3. 
In the clustered case Eq.p must be corrected by taking into account the 
fact that the softest gluon, in spite of the fact that is emitted in H, can be 
nevertheless clustered with the hard jets, and therefore does not contribute 
to Eo.p. Even for primary emissions this constitutes a LL contribution. 
The resummation of these multiple emission effects can be performed with 
the same MC used for non-global logs, and is the last bit that is needed to 
resum the dijet rate with the kt algorithm at NLL accuracy. 
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Figure 3. The coefficient of (as/27r)^C^ for the difference between resummed and NLO 
predictions for <j~^daldL, as explained in the text. Jets are clustered with the kt 
algorithm with /? = 1, and H is a rapidity slice with A 77 = 1. C 2 is the correction that 
needs to be applied at order to take into account multiple emission effects. 


4. Phenomenology of dijet observables 

Instead of going through the details of the resummed calculation for the 
dijet rate in DIS, we discuss what information on QCD dynamics we can 
gain from allowing symmetric Et cuts. Since in the region A —> 0 we are 
almost in a three-jet configuration, this measurement is complementary to 
three-jet event shapes, and allows one to investigate the coherence prop¬ 
erties of QCD emission from a three parton antenna. This observable has 
also the advantage that NP corrections are smaller than in event shapes 
distributions. This is mainly due to the fact that setting A —^ 0 does not 
put a direct veto on emitted parton transverse momenta, but rather on 
their vector sum. We have then emissions with large transverse momenta 
contributing also in the extreme A —> 0 region. This fact also affects the 
magnitude of non-global logs, which here^^ give a correction of order 10%, 
while for event shapes^ their contribution can be as large as 30%. 

One can consider other observables that have the same resummation as 
the dijet rate and can be more easily handled theoretically or experimen¬ 
tally. Among these are the distributions in the transverse energy difference 
A = Eti — Et 2 and in the azimuthal angle A(/> between the jets. For any 
of these variables V, one studies 


1 da 
a dV 



( 9 ) 
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Figure 4. Preliminary plots for the transverse energy differential distribution in DIS. 
The parameter X is used to estimate NNLL effects. 


As already stated, no further effort is needed to resum these distribution, 
since S(F) = <7q{Q)W{V), where W is the same as in eq. (4). Since from the 
previous analysis we have seen that S(V) ~ V for small V, one expects any 
of these differential distributions to approach a constant for V —> 0. This 
is what is already seen in the data for the A;/) distribution in hadronic dijet 
production at the Tevatron^®. This observable is particularly interesting 
because it gives access to the QCD fifth form factor, which is present when 
there are at least four hard emitting partons^^. 

A plateau for small A is also seen in the differential distribution in the 
dijet transverse energy difference in DIS computed with HERWIG^®. Fig. 4 
shows the comparison between the resummed result and that obtained from 
HERWIG. In this case jets are identified with the kt algorithm, and in 
the resummation non-global logs have been approximated with their value 
for the cone algorithm, which has been taken as an upper bound. The 
pretty good agreement between the two predictions, together with the fact 
HERWIG does not fully contain the matrix elements giving rise to non- 
global logs, suggests that the observable is dominated by soft and collinear 
emissions. 

Another interesting possibility could be that of making Eti—Et 2 global 
by modifying the particle recombination scheme. Actually, in a massless 
Eq scheme, one finds that emissions clustered with the jets do contribute 
to Eti — Et 2 , so that one can use the program GAESAR® to obtain auto¬ 
matically the resummed expression for its distribution. 
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The concluding remark of this overview is that pushing measurements in 
phase space regions where fixed order predictions are not expected to give 
an accurate description of the data opens up the possibility to investigate 
properties of QCD otherwise unaccessible. In particular, cross sections for 
jets with almost equal represent a yet poorly explored field, which 

could be complementary to the traditional event-shape measurements. 
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